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We study the entanglement distillability properties of thermal states of many-body systems. Fol- 
lowing the ideas presented in [D. Cavalcanti et al.. larXiv:07 05.3762 . we first discuss the appearance 
of bound entanglement in those systems satisfying an entanglement area law. Then, we extend 
these results to other topologies, not necessarily satisfying an entanglement area law. We also study 
whether bound entanglement survives in the macroscopic limit of an infinite number of particles. 
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I. INTRODUCTION 

Quantum Information Theory was born as a frame- 
work capable of describing the physics behind the pro- 
cessing, storage, and exchange of information under the 
rules of Quantum Mechanics [l|. Since it uses a very 
general and rather abstract approach, we do not need to 
talk about specific systems, interactions, or Hamiltoni- 
ans, but instead we can play with qubits, channels, logic 
gates, etc. Beyond information purposes, this abstraction 
has also been proven very useful when studying relevant 
questions in other subareas of physics, such as Condensed 
Matter Q , Statistical Mechanics ^] , Quantum Optics [1] , 
or Astrophysics 

In a recent work we have discussed the entangle- 
ment distillability properties of thermal states of some 
quantum many-body models with local interactions. We 
have shown the existence of a temperature range for 
which no pure-state entanglement can be distilled from 
the system despite being entangled. This type of irre- 
versible quantum correlations is also known as bound 
entanglement This result, which can be valid for 
systems of arbitrary size, was connected to the so called 
entanglement- area laws, a typical feature of these sys- 
tems that relates the entanglement of two distinguished 
regions to the area between them In the present 
contribution, we extend these ideas by considering two 
ways of addressing area laws and discussing the role they 
play in the appearance of bound entanglement. We also 
show that bound entanglement can appear in systems 
with different topologies, not necessarily fulfilling an en- 
tanglement area law. 

Let us start by saying few words about the main sub- 
jects considered here: distillability and area laws. It is 
well known that some entangled mixed states have the 
property of being distillable. This means that if one con- 
siders many copies of such a state, it is possible to purify 
the entanglement into a smaller number of entangled pure 
states using local operations and classical communication 
(LOCC) Although all two-qubit and one-qubit- 

one-qutrit states are distillable, there are states in higher 
dimension for which no LOCC operation can purify en- 
tanglement 01 ■ These states are called bound entangled. 
Specifically, in the case of a multipartite system an en- 
tangled state of n parties is bound entangled whenever 
the n parties cannot distill any pure-state entanglement 



out of it by LOCC. 

The first criterion able to detect bound entanglement 
was given by the Peres criterion : all distillable states 
have a non-positive partial transposition. Thus, find- 
ing an entangled state with positive partial transposition 
(PPT) guarantees non-distillability. In this paper we will 
use a quantitative version of this criterion, namely the 
negativity E]^{p) [l^l, to analyze the distillability prop- 
erties of thermal states. The negativity is defined as fol- 
lows: 
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Ai<0 



(1) 



being Ai the eigenvalues of p^^ , the partial transposition 
of p with respect to a given part A of the system. In other 
words, E]\r{p) is given by the sum of the absolute values 
of the negative eigenvalues of p^^ . So if the negativity of 
a state is zero, its partial transposition is positive, and 
the system is non-distillable (either separable or bound 
entangled), i.e., 

En{p) = p is non-distillable, 



En{p) > ^ p is entangled. 

Another useful and related quantity used throughout 
this paper is the logarithmic negativity Ei{p), given by 
Ei{p) = log2(l + En{p)). Clearly the implications above 
also apply to Ei [p) . 

Consider a system in a pure state and a bipartition 
of it into two complementary subsystems. It is well 
known that for bipartite pure states the entropy of one of 
the reduced systems uniquely measures its entanglement 
[1 ■ [T3|. As in Thermodynamics, one may expect that the 
entropy increases with the volume of the reduced state, 
but curiously this is not the case for the ground state of 
many models considered so far 0, [Hill, O Hill, 
In fact, it is instead seen that the entanglement be- 
tween the two complementary regions scales at most as 
their boundary area (in non-critical situations), a behav- 
ior that is generally named entanglement-area law. The 
above mentioned works approach entanglement-area laws 
in a variety of ways. Here, we will in particular focus on 
two different approaches. On one hand, we will address 
the entanglement-area relation by keeping fixed the size 



2 



of the system while changing the geometry of the biparti- 
tion. For example, in the case in which the distinguished 
subsystem is contiguous and the total system is consid- 
ered in the macroscopic limit such an approach is usually 
named as block entropy. When a system obeys an area 
law under this approach we will say that an area law of 
type I is fulfilled. On the other hand, one can consider 
also the opposite approach. Namely, how the entangle- 
ment scales when a given partition is kept fixed while 
changing the size of the system. An example of such 
an approach is given, again for contiguous partitions, in 
Ref. [l9|, where an area law in the half-half partition (i.e., 
when a contiguous group composed by half of the par- 
ticles belongs to the first subsystem and the other half 
to the second one) is established. In this case we will 
say that an area law of type II is obeyed. We consider 
both alternatives here, showing in which sense they are 
not equivalent and how they help in enlightening the ap- 
pearance of bound entanglement. 

The paper is organized as follows. In Sec. II we will 
point out the role that area laws play in the appearance 
of bound entanglement at finite temperatures. Then, in 
Sec. Ill we will focus on harmonic oscillator systems and 
analyze the emergence of different forms of area laws by 
studying different partitions and system sizes. We also 
consider a configuration where an area law is not seen in 
its simpler form and discuss the existence of bound entan- 
glement for this case. In Sec. IV we extend our analysis 
to the emergence of bound entanglement also for spin— i 
systems, obtaining very similar results as for harmonic 
systems. Sec. V is devoted to concluding remarks. 



II. BOUND ENTANGLEMENT AND AREA 
LAWS 

In Ref. [s] we considered systems exhibiting area laws 
and suggested that they are good candidates for pre- 
senting thermal bound entanglement in the macroscopic 
limit. This comes from the fact that when we increase 
the systems' size the ground-state negativity for some 
partitions increases (e.g. for the even-odd cut, where 
even particles belong to one subsystem and odd parti- 
cles to the complementary subsystem - see [13] ) while for 
other type of partitions it saturates (e.g. in the half-half 
geometry). It is then natural to expect that, when tem- 
perature is added to the system, the negativity in the 
half-half partition vanishes for a lower temperature than 
in the even-odd geometry. In other words, one expects 
that different partitions of the system become PPT at 
different temperatures (that we call threshold tempera- 
tures). An important feature of the considered systems is 
that they are translationally invariant and then all half- 
half partitions are equivalent. This observation implies 
that when the negativities of the half-half partitions are 
null no pure-state entanglement can be distilled by LOCC 
Q. On the other hand, there is still a temperature range 
where the negativity for the even-odd partition is strictly 



positive, which is enough to prove bound entanglement. 

As already mentioned, a lot of efforts have been de- 
voted to the study of entanglement-area laws for the 
ground state of various physical systems. Much less is 
known about equivalent laws for the case of finite tem- 
perature, where the state of the system is in a thermal 
mixture. Partly, this is due to the complexity character- 
izing the structure of entanglement for mixed states. For 
example, one of the known exact results concerns again 
the entropy of a contiguous subsystem [2l|. However, for 
mixed states, this quantity is no longer a measure of en- 
tanglement. To the best of our knowledge, the most gen- 
eral result at finite temperature has been derived recently 
by Wolf et al. [23| and gives a bound to the mutual in- 
formation between two complementary subsystems. We 
recall that the mutual information is a measure of the to- 
tal amount of correlations, both classical and quantum. 
Hence, it trivially gives an upper bound to the entan- 
glement. An important feature of such a bound is that 
the dependence with the temperature T and the area A 
between the complementary regions is factorized (in par- 
ticular the bound scales linearly with A and T^^). On 
the basis of this result, one may argue that a trivial area 
law, i.e., in a factorized form, holds for the negativity 
Eiq whenever: 

En < I{T)9{A) , (2) 

where /(T) and g{A) are generic functions also depend- 
ing on the parameters of the system. In particular, we 
assume that these functions do not depend on the way 
of partitioning the system. As said, this type of rela- 
tion holds for the bound found in Ref. [l^l, in the case 
of the mutual information. In case the inequality in the 
previous formula becomes an equality, we refer to this 
form of area law as a strict area law. Notice that such a 
form of area law has been shown to hold for a nearest- 
neighbor harmonic ring at finite temperature, in the case 
of even-odd partition [6|. There, numerical evidences of 
the validity of a strict area law have been reported also 
for half-half partitions as well as for the analogous cases 
in a spin ring. 

Let us now consider how a strict area law affects the 
existence of bound entanglement. As said, the key ingre- 
dient in the recipe above in order to show the presence 
of bound entanglement is that different partitions of a 
systems become PPT at different temperatures. As a 
consequence, if a strict area law of type I holds then no 
bound entanglement should be expected, since all par- 
titions become PPT at the same temperature, namely 
when /(T) = Notably, we have found no sys- 

tem showing this behavior, even considering models that 
exhibit the same entanglement for different partitions in 
the ground state, as we will report in detail in the next 
Sections. On the contrary we observed, for any system 
taken into account, that a strict area law of type I does 
not hold and that different partitions become PPT at 
different threshold temperatures (in any case, the bound 
given in Ref. [l^] is of course not violated). This means 
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that for a fixed system size, there is a temperature range 
for which bound entanglement is present. 

Now, one may wonder whether this holds also in the 
macroscopic limit of an infinite number of particles. In 
this case, it is the validity of a strict area law of type 
II that does give the key to a positive answer. In fact, 
it ensures that the threshold temperatures for each of 
the partitions chosen to reveal the presence of bound en- 
tanglement {e.g., even-odd and half-half partitions) stay 
constant as the size of the system increases. In particu- 
lar, it ensures that the range of temperatures for which 
bound entanglement is present survives up to a macro- 
scopic level. We will see explicitly that this is actually 
the case for some models consisting of harmonic chains 
and that the same behavior seems to be valid also for 
spin chains. In order to show that different scenarios 
may arise in other topologies, we also considered systems 
with a star configuration. There, we will see that a strict 
area law of type II is no more valid. The next sections 
are devoted to show explicitly the ideas explained above. 

III. HARMONIC OSCILLATORS 

Consider a system consisting of N harmonic oscillators 
interacting via the following Hamiltonian: 

where Xi and "Pi represent position and momentum oper- 
ators for each oscillator respectively (i = 1, . . . , N). 

The matrix V describes both the on-site interaction 
(given by the diagonal elements) and the coupling be- 
tween oscillators i and j (non-diagonal terms). This 
Hamiltonian is quadratic in the canonical coordinates 
and the oscillators are coupled through their position 
degrees of freedom which sets both the ground and 
the thermal states to be Gaussian. In this way these 
states are completely determined by their covariance 
matrix 7 defined as follows. Take the vector S — 
(xi, . . . , xn,Pi, . . . ,Pn), we then have 

7fc, = Re {TT{Q[Sk - Sk] [Si - Si]}) (4) 

where g is the density matrix of the state and Sk — 
Tr(pS'i;). If we consider the thermal state g = 
exp[-iJ/r]/Tr{exp[-i7/T]} at temperature T, the cor- 
responding covariance matrix is given by p^j : 

7(T) = [V-^/^W{T)] e [V^''^W{T)] , (5) 

where 

T/F(r) = ljv + 2[exp(t/i/VT)-lAr]-\ (6) 

and Iat denotes the NxN identity matrix. In the ground 
state case 1^(0) is given by the identity matrix and so 

7(0) = -1/-1/2 e 1/1/2. 




FIG. 1: Pictorial representation of the considered configu- 
rations. A. Particles interacting via nearest-neighbors cou- 
plings. B. Star configuration: a central particle interacts 
equally with the rest of the particles. 

An analytical expression for the entanglement (quan- 
tified by the log-negativity Ei [l2|) between two com- 
plementary groups of oscillators, A and B, was given in 
terms of the covariance matrix of the state, which can be 
written, in turn, only in terms of the matrix V. Then 
one gets the general formula for the log-negativity of a 
thermal state at temperature T: 

= ^log2{max[l,Afc(Q)]}, (7) 

k=Q 

where Q = P uo- Plo+ and w± = W{T)-W^^. We 
denoted by {A/c[Q]}^Jo^ spectrum of the matrix Q, 
whereas P is aii N x N diagonal matrix with the i-th 
entry given by 1 or —1 depending on which group, A or 
B, oscillator i belongs to. 

In Ref. [sl we considered the thermal states of Hamil- 
tonian ([3|) with a circulant potential matrix V given by 

y = circ(l,-c,0,...,0,-c). (8) 

i.e., the particles interact via nearest-neighbors interac- 
tions (see Fig. [IK)- We analyzed the entanglement be- 
havior of the even-odd and half-half partitions while the 
system's size N is increased. A strict area law of type 
II for the log-negativity was analytically obtained for 
large N, specifically the entanglement for a given par- 
tition changes proportionally to the area A . The change 
of the system temperature just affects the rate the en- 
tanglement increases with iV for an even-odd partition 
(^ = in this case) and the entanglement saturation 
value for the half- half partition (^ = 2). In particular, 
the threshold temperatures does not depend on A , com- 
patibly with a strict area law of type II. 

We now investigate the area law of type I, that is for 
fixed N — Ntot and varying the area by changing the 
geometry of the partition. Let us focus on two different 
ways of partitioning the chain. First, we consider the par- 
tition of the system consisting of Ntot = 2" particles into 
2"' alternate blocks. The area associated to such a par- 
tition is given by ^ = 2"''. For rih — n one retrieves the 
even-odd partition and for Uf, — I the half-half partition. 
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FIG. 2: Log-negativity as a function of the area in the case 
of nearest-neighbor harmonic ring of fixed size Ntot = 2^ and 
c = 0.4. From top to bottom the inverse temperature /3 = 1/T 
is given hy [3 = 2.5, 2.4, 2. The system has been partitioned 
into symmetric alternate blocks (see text for details). 
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FIG. 3: Same as Fig. O but partitioning the system in a non 
symmetric way (see text for details). The inverse tempera- 
ture 13 = 1/T is given by /3 = 1.87, 1.865, 1.863, from top to 
bottom. 



In Fig. [5] we depicted the log-negativity as a function of 
the area for Ntot = ^ and c = 0.4 at different tempera- 
tures. One sees that, for some fixed temperatures, there 
are partitions with PPT while others with non positive 
partial transposition, meaning that not all of them be- 
come PPT at the same temperature. A strict area law 
of type I is then violated, allowing for the presence of 
bound entanglement. 

We also considered another way of gradually partition- 
ing the system, again for fixed system size. Starting from 
an even-odd partition i^A = Ntot) we took one particle at 
a time, and transferred it from, say, the even block to the 
odd block. The process ends with one of the blocks com- 
posed by only one particle (^ = 2). In this way the area 
is decreased by 2 at each step. In Fig. [3] we depicted the 
log-negativity as a function of the area for Ntot = 100, 
c = 0.4 and different temperatures. Again one sees that 
by increasing the temperature, some partitions become 
PPT while others are not. 

Through these examples we can see that a strict area 
law of type I does not hold. However, notice that even 
if the entanglement does not vary strictly proportionally 



to the area of the subsystem, it still increases with the 
area. The details of such a behavior of course depend on 
many factors, such as the way the system is partitioned 
and the entanglement measure. Nevertheless a general 
feature seems to be independent of the partitioning: the 
more the interaction bonds intercepted by the partition 
the more the entanglement across it. 

Summarizing, we have seen that in the case of a har- 
monic nearest-neighbor ring a strict area law of type I 
is violated, allowing for the presence of bound entangle- 
ment, whereas a strict area law of type II is valid, ensur- 
ing that bound entanglement survives for large systems. 
In the next section, we show that a similar behavior is 
also found for spin rings. 

Let us now turn to a system in which a strict area law 
of type II is not valid either. As anticipated, we studied a 
different topology, namely a star configuration (see Fig[T] 
B). The system is described by the Hamiltonian ^ with 
potential matrix given by ^1^1 = 1-|-(7V— l)c, F^f — 1 + c, 
V-l^ = -c and Vif = otherwise {2 < i < N , c> 0), i.e., 
all the oscillators are equally connected to a central one. 
Clearly, translational symmetry does not hold anymore. 
The area law of type I is violated also in this case, al- 
lowing to find a temperature range in which the state is 
bound entangled. In Fig.|3]we depicted how the threshold 
temperatures for which the log-negativity is zero, Tj^i^^ 
(half-half partition) and T^f^" (central particle versus the 
outer ones) , vary with N . In the region between these two 
curves, bound entangled states are present. Nonetheless 
this is not enough to guarantee that bound entanglement 
survives for a large number of particles. First, notice that 
the range of temper ature Tj^il ^ - Tg° is no longer constant 
with iV, namely, a strict area law of type II is not valid 
in this case. Second, our numerical calculations suggest 
that the entanglement between the central particle and 
the rest (which is the largest for this configuration) goes 
to zero as TV ^ oo. Actually, an analytical expression 
for the log- negativity can be guessed for T = 0. In this 
case, being the state pure, the information about the en- 
tanglement is completely given by the reduced covariance 
matrix 7i.ed of the central particle. The latter is simply 
given by the elements of 7(0) = y~i/20yi/2 correspond- 
ing to the central particle itself. By calculating explicitly 
7red for a small number of particles A'^, one can recognize 
the following structure: 

/I N -I r, — \ f I N -I \ 

(9) 

The negativity Eiq between the central particle and the 
rest is now simply given by 

En = max[0, ^-^] , (10) 

where v — ^fK — \/A — 1 and A is the determinant of 
7rod (see, e.g., Ref. [23]). Assuming now that the struc- 
ture given in Eq. ^ holds for a generic A'^ we can extrap- 
olate the behavior of the negativity in the macroscopic 



FIG. 4: T/^;° (solid line) and Tl^i'' (dashed line) as a function 
of the number of oscillators for the star system with in- 
teraction (see text for details). The coupling constant is 
c = 1 and similar behaviors have been found for different val- 
ues of c. We see that the range of temperatures for which the 
state is bound entangled varies with A'^ for this system. Inset: 
Log-negativity of the central vs. rest partition as a function 
of 1/N, for c — 1 and /3 = 1. Notice that for large A'' the 
log- negativity tends to vanish (see text). The discontinuity 
simply comes from the fact that A'' is integer. 

limit. In fact, since A ^ 1 in the limit N ^ oo, the neg- 
ativity itself goes to zero. Considering now the generic 
case at T > 0, it is then reasonable to expect that the 
log-negativity goes to zero too for a large number of par- 
ticles. As said, our numerical calculations confirms this 
intuition (see also the inset in Fig.S]). As a consequence, 
the system is fully PPT in the macroscopic limit. Such a 
configuration gives then a non-trivial example for which 
the absence of a strict area law of type II does not al- 
low to clearly identify a range of temperatures for which 
bound entanglement is present in the macroscopic limit 

IV. SPIN SYSTEMS 

The scope of this section is to extend the previous anal- 
ysis to spin systems. We concentrate on the thermal state 
of systems composed by N spin-^ particles, interacting 
with the Hamiltonian 

N 

H = - + + (11) 

<i,j> i=l 

where the pairs of indices i and j over which we sum 
define the topology of the system. 

In the nearest-neighbor configuration (see FiglTJ^) we 
proceeded as we did for FigO and progressively increase 
the boundary area between two regions, starting in the 
half-half partition and changing particles from one par- 
tition to the other. In this situation, for a chain with 
N = 10, we can observe (FigE]) again the presence of 
bound entanglement, since we have that some partitions 



FIG. 5: For a system of A'^ = 10 spins with nearest-neighbor 
interaction with hamiltonian given by Eq.((TT} and h — 1.9, 
we plot the negativity vs. the area of the boundary between 
the two regions, for T = 3, 3.15, 3.25. For certain values of T 
we only observe PPT entanglement above a certain number 
of bonds connecting the two regions. 

are entangled while others are PPT at the same temper- 
ature. This again means that a strict area law of type I 
is not valid. On the other hand the validity of an area 
law of type II has been numerically shown up to twelve 
particles in Ref. Q , in analogy to what we have seen for 
the harmonic ring. This features strongly support the ex- 
istence of bound entanglement in the macroscopic limit 
also for spin rings. 

Let us now move to a system with a star configuration 
as in FiglTjB, and consider the negativity corresponding 
to partially transpose either the middle particle or one of 
the outer particles. A remarkable feature of this system 
is that the ground state entanglement for both partitions 
is the same for any fixed N, actually these partitions are 
both maximally entangled. This can be easily seen by re- 
callirig the explicit expression of the ground state given in 
Ref. |24| . One may then wonder whether such a behavior 
holds also at non-zero temperature, which would suggest 
that a strict area law of type I is valid. Our calculations 
show that this is not the case. Again, different partitions 
become PPT at different temperatures, as can be seen in 
Figs. [H] and [T] Notice that the central particle now be- 
comes PPT at lower temperatures with respect to the ex- 
ternal ones. Nevertheless, one can conclude the presence 
of bound entanglement also in this case by computing 
the threshold temperature for a half-half partition. 

Another interesting feature of this system is that, by 
symmetry reasons, the entanglement between the middle 
particle and the external spins is independent of the sys- 
tem size (see Fig 15]). This is due to the fact that all the 
eigenstates of this system are of the form [23| : 

-^(|0)|a™,,)±|l)|a:„,^.)) (12) 

where the first ket correspond to the central particle, 
whereas the second one to the external particles. The 
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cles due to computational hardness. Although not shown 
here, similar results can be found for other types of in- 
teractions, e.g. using Heisenberg-type hamiltonians. 



FIG. 6: Negativity for a partition consisting of the central 
particle vs. the others as a function of the temperature T for 
a star system of A'^ = 4, 6, 8, 10 particles. The Hamiltonian is 
given by Eq. (|lip with /i = 0. We can see that the negativity 
for different A^'s coincides, showing its independence on the 
system size. 



states \am,j) and are orthonormal eigenvectors of 

a high dimensional fictitious spin. The key point here is 
that the partial transpositions with respect to the cen- 
tral particle do not change the structure of these eigen- 
vectors and this is true for any N . As a consequence, 
once expanded the thermal state in the eigenbasis above, 
one can see that T^f^ does not depend on TV. In other 
words a strict area law of type II holds for this partition. 
Notice however that this behavior does not hold in the 
case of other partitions, as we can see in FiglT] for the 
case of one external particle with respect to the rest. In 
particular, the threshold temperature for this partition 
increases with the system size. This fact, considering 
that the threshold temperature is size-independent for 
the middle particle partition, suggests that the tempera- 
ture gap for which bound entanglement appears between 
the two partitions increases with the system size. Recall 
that this gap appeared constant between the half-half and 
even-odd partitions @ in systems with nearest-neighbor 
interaction. 

We can also check the negativity corresponding to 
transpose one of the external particles, for different sys- 
tems sizes at a fixed temperature. In this scenario, one 
can see a peculiar behavior: for low temperatures the 
entanglement for the central vs. the rest partition can 
decrease as the system increases, but for higher temper- 
ature the opposite holds. The crossing temperature ap- 
pears clearly non-trivial, between 2.2 and 2.4, see the 
inset of Fig[71 Again the entanglement in different par- 
titions do not vanish at the same temperature, allowing 
the presence of bound entanglement. 

Before ending this section, let us stress that our numer- 
ical calculations are restricted to small number of parti- 




FIG. 7: The negativity for a partition containing only one 
of the external particles vs. the others as a function of the 
temperature T for a star system of A'^ = 4, 6, 8, 10 particles, 
where the interactions are described by Ea. (lll[l with h = G. 
The inset shows in detail the crossing section, and the different 
temperatures where the partition becomes PPT. 



V. CONCLUSIONS 

To conclude, we have extended the results of Ref. @ 
and considered different ways of studying the entangle- 
ment distillability properties of thermal states of many- 
body systems. We have considered systems of harmonic 
oscillators and spin-one-half particles in a chain and star 
topology. In general, our results show that a strict entan- 
glement area law of type I (when changing the partitions 
for a fixed system size) is not fulfilled. Since the dif- 
ferent partitions become PPT at different temperatures, 
bound entanglement appears for a temperature range in 
a natural way. Concerning the preservation of this range 
of temperatures when the system size is increased, we 
pointed out that a different approach to area laws should 
be addressed. In particular, an entanglement area law 
of type II (when changing the system size) is then use- 
ful to prove the presence of bound entanglement in the 
macroscopic limit of an infinite number of particles. 
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